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Abstract 

By combining the definition of the Wigner distribution function 
(WDF) and the matrix method of optical system modeling, we can 
evaluate the transformation of the former in centered systems with 
great complexity. The effect of stops and lens diameter are also con- 
sidered and are shown to be responsible for non-linear clipping of the 
resulting WDF in the case of coherent illumination and non-linear 
modulation of the WDF when the illumination is incoherent. As an 
example, the study of a single lens imaging systems illustrates the 
applicability of the method. 

1 Introduction 

The Wigner distribution function (WDF) provides a convenient way to de- 
scribe an optical signal in space and spatial frequency [TJ Ej . The prop- 
agation of an optical signal through first-order optical systems is well de- 
scribed by the WDF transformations |S1 IH E] , allowing the reconstruction 
of the propagated signal. Real optical systems are not first order and the 
use of the WDF for optical system design presumes the ability to predict 
how it is transformed by systems with aberrations. Almeida jH] proposed a 
method to determine the aberration coefficients for optical systems using ma- 
trix methods and calculated the necessary coefficients for 7th-order modeling 
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of centered systems based on spherical surfaces The extension of the matrix 
method to cylindrical surfaces has also been proposed [9 . Based on the fact 
that the WDF lies between Fourier and geometric optics, we show that geo- 
metric optics matrix coefficients can be used to predict WDF transformations 
and hence can play an important role in optical system design. 

2 Transformation of the WDF 

The Wigner distribution function (WDF) of a scalar, time harmonic, and 
coherent field distribution <£>(q, z) can be defined at any arbitrary z = const, 
plane in terms of either the field distribution or its Fourier transform Tp(p) = 



where q is the position vector, p the vector of the conjugate momenta, k = 
2tt/\ and * indicates complex conjugate. In the present work we will be using 
mainly quasi-homogeneous light, in which case the WDF can be defined as 



where i(q) is a non- negative function which we call the intensity and s(p) 
is the Fourier transform of the positional power spectrum s(q) and is also 
non-negative. 

If the position coordinates are x, y, z and the ray direction cosines are 
u,v,w, the position and conjugate momenta vectors are given by 



/ y?(q) exp(-ikq r p)dq 0H2]: 
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where n is the refractive index of the optical medium. 
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We will assume that the optical system is characterized by a transfer map 
between the initial phase space coordinates, q J , p* and the final ones, q^, p^. 
If M. represents the transfer map: 



q' 




(6) 



We will also write expressions like q^ = M.a£ or x* = M.x l to represent the 
dependencies of each of the final coordinates on the original ones. 

The transfer map can always be inverted; a simple physical argument is 
sufficient to prove it: The transfer map is the relationship between the ray 
coordinates on the input plane (p*, q l ) and the corresponding coordinates 
on the output plane (p^,q^); if two rays share the same coordinates on the 
output plane they are the same ray and so it is always possible to map the 
output onto the input. This type of reasoning is valid within the scope of 
geometrical optics, which corresponds to the conditions for existence of a 
transfer map. 

Having established that the map can be inverted the WDF transformation 
is governed by the equation 

W f (q f ,p f ) = |/i(q / ,p / )|W i (A^- 1 q / ,>f- 1 p / ), (7) 

where the factor \h(qf ,pf)\ accounts for the energy conservation between 
input and output and is the ratio between the elementary hypervolume in 
input phase space and the corresponding mapped hypervolume in output 
phase space: 



|Zi(q',p^ 




2 dx l dy t du t dv t 



dx-f dyf duf dvf 

if J is the jacobian of the map transformation we can write [Zj 



(8) 



\h{q?,pf)\ = ± (9) 

Eq. (J2J) is of special interest when the transfer map can be expressed in 
closed form, which is the case if matrices are used jHHOHOI- I n this situation 
the output coordinates are expressed as polynomials in the input coordinates 
or vice-versa. Almeida |B] showed that this method can be extended to any 
desired degree of approximation, at least for centered systems, and published 
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the coefficients for the 7th-order matrices of systems based on spherical sur- 
faces. 

There are two methods of map inversion in matrix optics. The first one is 
a straightforward matrix inversion and can be used in many cases; the second 
one, applicable in all circumstances, consists in reversing the optical system 
and recalculating all the matrix coefficients. We can thus find Ai, and 
J for any centered optical system, no matter how complex. There remains a 
question about the aperture stops which is dealt with below. 

In order to model a system with matrices we start by defining a gener- 
alized ray of complex coordinates Q l = x l + jy l and P l = n(u l + jv l ); this 
ray is described by the 40-element monomials vector Q l &, built according 
to the rules explained by Kondo et al. jH] and Almeida [U]. If the ray is 
subjected to a transformation described by matrix M, then the output ray 
has coordinates (Qf,Pf) and is represented by the monomials vector Q^&, 
such that: 

Q / & = MQ i &. (10) 

For an axis symmetric optical system, in the 7th-order, matrix M will 
result from a product of 40 x 40 square matrices with real elements. Each 
matrix in the product describes a specific ray transformation. The elementary 
transformations can be classified into four different categories: 

Translation: A straight ray path. 

Surface refraction: Change in ray orientation governed by Snell's law. 

Forward offset: Ray path between the surface vertex plane and the sur- 
face. 

Reverse offset: Ray path from the surface back to the vertex plane, along 
the refracted ray direction. 

The ray itself is described by a 40-element vector comprising the mono- 
mials of the complex position and conjugate momenta coordinates that have 
non-zero coefficients; the first two elements of this vector are just the complex 
coordinates (Q,P). 
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Considering Eq. (JTUJ) the mapping equation (JUJ) takes the form 



Q f = V Q (Q\Q l *,P\Pn; 

P f = V P (Q i ,Q i *,P i ,P i *); 

Q l = VQ\Q f ,Q f *,P f ,P f *); 

P i = Vp\Q f ,Q f \P f } P f *). 



(11) 
(12) 

(13) 

(14) 



The symbols Vq, Vp, Vq , Vp mean polynomial expressions of the variables 
in parenthesis. Eqs. (fE^ and (fTl| can now be used to evaluate Eq. first 
and then Eq. (J7|). 

3 Stops and pupils 

Any system analysis is incomplete without consideration of the effect of the 
various stops along the optical path; this analysis cannot be incorporated 
in the matrix description and deserves special treatment. Paraxial theory 
tells us [TU] that we can find one most limiting stop whose images in object 
and image space are known by entrance and exit pupils, respectively. The 
theory goes that the entrance pupil establishes the width of the beam entering 
the system while the field angle is established by the second most limiting 
stop imaged onto object space; the images of the same stops in image space 
set corresponding limits to the rays leaving the system. It is not necessary 
to leave paraxial theory to find that these concepts are insufficient for the 
complete description of the beam constraints within the system and we are 
led to the concept of vignetting. 

Moving from the paraxial approximation to high-order the problem in- 
creases in complexity and even the concepts of entrance and exit pupil loose 
significance in view of the high aberrations present when an internal stop is 
imaged to either object or image space [TT]. Ray tracing software usually 
avoids the problem by imposing restrictions as the rays cross each stop's 
plane [T2] . 

In order to tackle the problem in phase space, we will define scene as an 
optical field distribution that spans — oo < |q| < +oo in space coordinates 
and — n < |p| < n in conjugate momenta. A scene cannot contain com- 
ponents with |p| > n because these components originate evanescent waves 
that are considered faded out [131 E] . What the paraxial theory says is that 
the entrance pupil clips the scene in p coordinates, while the second most 



5 



limiting stop is responsible for clipping in q coordinates. If only the merid- 
ional plane is considered, to reduce the dimensions to 2, the stops create 
a parallelogram area in phase space, with two sides parallel to the p axis, 
where the WDF is non-zero. Vignetting must be seen as a departure from 
that form, meaning that, for the extreme values of q, the angular spread of 
the rays may be different from the central one. 

In order to understand the stop effects on the WDF we consider it as a 
modulator, in which case the following relation applies [To] : 

W"(q, P) = ^ / W m (q, p - pO W*(q, p'W, (15) 

where W m (q, p) is the WDF of the modulating function m(q). Eq. (JTSJl rep- 
resents a two-dimensional convolution of the Wigner distribution functions 
W m (q,p) and W™(q, p) with respect to the frequency variables and a mere 
multiplication with respect to the space variables. 

A stop is a special kind of modulator. In coherent illumination the stop 
has a modulating function that equals unity within the stop area and is zero 
elsewhere. Furthermore, as we are usually dealing with stops that are very 
large compared to the wavelength, Eq. (fT5|) results in clipping of the local 
WDF in the space domain: 

W f (q, p) = W l (q, p), within the stop, 

W / ^( ( l?p) = 0, elsewhere. (16) 

In incoherent illumination the stop modulating function is the auto-correlation 
function of the stop transmittance function [TOl IT3] . If as before the stop di- 
mensions are large compared to the wavelength, Eq. (|15|) can be written 

W f (q, p) = W l (q, p)S(q), within the stop, 

W(q>p) = 0> elsewhere. (17) 

The stop auto-correlation function is defined by 

S(q)= J 5(q + q')5(q)dq'. (18) 

The translation of the stop modulation onto an equivalent effect of the 
original scene's WDF depends on the sort of transformations the latter has 
incurred up to that point. When the signal encounters the first stop in 
the system the only transformation that the WDF has suffered is a spatial 
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shearing, which is linear in the paraxial approximation and non-linear if 
wide angles are considered ^3]. If the distance from the scene to the stop 
is large the angle subtended by the stop will be virtually independent from 
the position coordinates on the scene and the stop effect will be virtually 
equivalent to a clipping or modulation on the spatial frequency domain. This 
is what an entrance pupil is supposed to do and so we state that an entrance 
pupil is a concept valid in the paraxial approximation, when the scene is very 
far from the optical system. 

The effect of further stops along the system is more difficult to under- 
stand. Let us assume that we are dealing with small angles, such that parax- 
ial approximation is indeed applicable, that we have converted the existing 
stops to their equivalents in object space and let us consider just the two most 
significant ones. The problem has been reduced to free-space propagation, 
characterized by linear shearing of the WDF. 

Fig. [T] illustrates the situation described above; object point Pi is an 
axis point and obviously stop S2 is the entrance pupil, responsible for the 
limitation on the rays that enter the system. According to general practice, 
we would define the field limits as the points on the rays that pass on the 
edges of the stop Si and the center of the entrance pupil; point P2 is one such 
point. Naturally, stop Si also introduces limitations on the rays that enter 
the system, besides its prime function as field limiter; this effect is known as 
vignetting. 

The effect in phase space is illustrated in Fig. |2] where a similar situation 
is depicted. Fig. |21 a) shows the effect of the entrance pupil on the local 
WDF with a clipping to its own width. The field pupil, stop Si, produces 
its clipping on the WDF back-propagated from the entrance pupil, as shown 
in Fig. 121 b), and the resulting double-clipped WDF is back-propagated to 
the scene plane, as shown in Fig. |2]c). The resulting parallelogram shape 
is the representation in phase space of the signal that can, in fact, enter 
the system; it is clear that for a point on the axis, q = 0, the stop S2 is 
responsible for determining the admittance angle, while stop Si is, to a great 
extent, responsible for determining the dominion of q, which is exactly what 
we call field. The vignetting effect is visible for extreme values of q, for which 
stop S2 no longer determines completely the admittance angle. 

The extension of the above procedures to a general mapping situation, 
outside the paraxial approximation, must be done carefully. In the next 
section we study a complete mapping situation, illustrating both the WDF 
transformations and stop consideration. 
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4 Example 



The case below was chosen not for its particular applicability but for its 
ability to demonstrate and highlight the possibilities opened by the matrix 
mapping and WDF used together. 

We will consider a simple imaging system composed of a single convex 
lens and a field stop on the image plane. The lens was chosen to produce a 
high degree of aberrations, so that the non-linear effects are clearly visible. 
The lens is plano-convex, with the flat surface facing the image plane, and 
has a refractive index of 1.56; the convex surface has a radius of 6.5 x 10~ 2 m 
and the central thickness is 2.0 x 10 -2 m. The lens and field stop diameters 
will be decided later on, upon examination of the aberrations present in the 
image. 

Eqs. (jlll to inj) were established for the system in consideration using the 
method outlined by Almeida [H] using Mathematica ^HI- The same software 
package was also used for all the further calculations. The analysis was 
carried out on a meridional plane, so the detection of aberration effects such 
as astigmatism is out of the question. The input scene was defined according 
to Eq. in terms of its WDF as: 

= 1 + sin-^, (19) 

where I is a parameter used to control the detail on the scene; all the graphics 
were plotted with I — 5 x 10~ 2 m. 

The input scene was located at 13.7 x 10~ 2 m, so that the image was 
formed at 5 x 10 -2 m. Fig. |3] shows the input scene in phase space and the 
output WDF. The input appears as a series of light and dark bands, showing 
the independence of the corresponding WDF on the p coordinate, character- 
istic of spatially incoherent light, a special case of quasi-homogeneous light 
[HE]. The output shows the same bands, reduced in width due to a magni- 
fication factor lower than unity and distorted by aberrations. A qualitative 
analysis of the aberrations is indeed interesting. 

The S-shape of the bands results from spherical aberration of various 
orders, with predominance of the third-order. The reduced width of the 
bands for higher values of q is characteristic of barrel distortion; this so high 
that the signal does not exist above \q\ > 0.5, except for the effect of spherical 
aberration. This is similar to a fish-eye objective. Field curvature is clearly 
visible as tilting of the central portion of the bands for high |g|. Coma results 
in an asymmetry of the S-shape. 
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Clearly we have performed a mapping with an infinite diameter lens, 
which only works mathematically. Considering the radius and width of the 
curved surface, we have established a lens diameter of 4 x 10~ 2 m. The 
lens diameter was given to a stop located on the vertex plane and the edges 
of this stop were mapped forward, through the lens and free-space, to the 
image plane, and backward to the input scene plane. The maps of the lens 
diameter stop are superimposed on the corresponding WDFs as dashed lines. 
We decided to use a field stop on the image plane, in order to limit the image 
to an area of low aberration; a diameter of 4 x 10~ 2 m was also chosen for 
this stop. The field stop was mapped onto the input scene plane and is 
shown as a solid line superimposed on both figures. If we were dealing with 
coherent light the area of both WDFs common to the zones defined by the 
two stops would be the area relevant for the image formation; in fact, to put 
it correctly, the image WDF should have been made equal to zero outside 
that area. 

For a one-dimensional stop Eq. (j!8|) becomes |X3) 



where d is the half-width of the stop. So, in incoherent illumination, rather 
than clipping the local WDF, the stop produces a gradual transition from 
full intensity to zero with twice the width of the stop. When propagated 
to either the image or the object planes this transition manifests itself as a 
gradual transition of the WDF from the full mapped value to zero guided 
but not delimited by the stops' traces, see Fig. HJ 

5 System analysis 

Although not presented in this paper, it would be possible to extract a lot of 
information about the system from the image WDF modulated by the stops. 
The field distribution would be obtainable directly from an integration of the 
image WDF in the variable p; the integration limits would be established by 
maps of stops twice the width of real the stops. It is clear that, within the 
region delimited by the field stop, there is a reasonable reproduction of the 
original scene. 



S(q) = 1 



S(q) = 




q\ < 2d, 
q\ > 2d., 



(20) 
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The point spread function for an input point go could be evaluated con- 
sidering a different input scene, such as W l (q,p) = S(q — qo), and again 
integrating the output WDF in p. The MTF could also be evaluated using 
the same scene but performing the integration in q. 

6 Conclusions 

The authors presented a method to evaluate the WDF transformations of 
an optical signal that passes through a system, in the context of geometrical 
optics. Using matrices it is possible to model centered systems up to any 
desired order of approximation; the authors have shown that the same matrix 
method can be used for the evaluation of the WDF transformations. 

The effect of stops and lens diameter could also be accounted for leading 
to the definition of clipping traces on both the input and output WDF and 
to the outline of methods to evaluate the resulting field distribution, point 
spread function and MTF. 
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Figure 1: The effect of stops on propagation angles. S2 is the entrance pupil 
and object point P2 is on the edge of the field. 
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Figure 2: Effect of stops in phase space, a) Clipping at the entrance pupil, 
b) the WDF is back-propagated from the entrance pupil to the field pupil 
and clipped, c) back-propagation to the scene plane. 
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Figure 3: Transformation of the Wigner distribution function through a lens; 
a) input distribution, b) output distribution. Both figures show the clipping 
effect of the lens diameter (dashed line) and the field stop (solid line). 
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Figure 4: The Wigner distribution function on the image plane shows the 
effect of stops. With incoherent illumination the stops don't clip the distri- 
bution but apply a smoothing from the center of the stop to twice the stop 
width. 
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